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Abstract 

Since the first TeV blazar Markarian (Mrk) 421 was detected in 1992, the number 
of established TeV 7-ray emitting BL Lac objects has grown to 6, with redshifts 
ranging from 0.031 (Mrk 421) to 0.129 (H 1426+428). The intensive study of these 
sources has had a major impact on our understanding of the blazar phenomenon. 
The most notable observational results have been extremely fast large amplitude 
flux and spectral variability on hour time scales, and a pronounced X-ray - TeV 7- 
ray flux correlation. In this paper we discuss recent observational results and report 
on progress in their theoretical interpretation. 
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1 Introduction — TeV Blazars 

The EGRET (Energetic Gamma Ray Experiment Telescope) detector on board 
of the Compton Gamma-Ray Observatory discovered strong MeV 7-ray emis- 
sion from 66 blazars, mainly from Flat Spectrum Radio Quasars and Uniden- 
tified Flat Spectrum Radio Sources (^2). Ground-based Cherenkov telescopes 
discovered TeV 7-ray emission from 6 blazars, 4 of which have not been de- 
tected by EGRET. All the established TeV detections belong to the class of 
BL Lac objects, blazars with relatively low luminosity but with Spectral En- 
ergy Distributions (SEDs) that peak at extremely high energies (see Table 1). 

The large detection area of Cherenkov telescopes of ~10^ m^ makes it possi- 
ble to sample the 7-ray lightcurves with a time resolution of several minutes. 
Large amplitude flux variability on 30 min time scales implies that the TeV 
emission originates from a small region very near the supermassive black hole 



Source 


z 


Discovery &: Confirmation 


Mrk 421 
Mrk 501 
lES 2344+514 
lES 1959+650 
PKS 2155-304 
H 1426+428 


0.031 
0.034 
0.044 
0.047 
0.116 
0.129 


Punch et al. 1992, Petry et al. 1996 
Quinn et al. 1996, Bradbury et al. 1997 
Catanese et al. 1998, Tluczykont et al. 2003 

Nishiyamaet al. 1999, Holder et al. 2003, Aharonianet al. 2003 

Chadwick et al. 1999, Hinton et al. 2003 
Horan et al. 2002, Aharonian et al. 2002a 



Table 1 

TeV blazars with a significant detection by at least 2 experiments (as of July 2003). 
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f|3^ : I47I ) . The relativistic Doppler Factor of the emitting plasma is defined as 
= r(l — (3 cos(^)), where F is the bulk Lorentz factor of the emitting 
plasma, /? is its bulk velocity in units of the speed of light, and 6 is the angle 
between jet axis and the line of sight as measured in the observer frame. The 
requirement that the emitting volume is optically thin for the TeV 7-radiation, 
together with assumptions about the optical to UV energy spectrum of co- 
spatially emitted synchrotron emission gives a lower limit of (5j > 9 for the 
emitting plasma fj3Qi). The Lorentz boost of the observed jet luminosity by the 
factor makes the rather weak jet emission readily detectable. 

The main scientific drivers for broadband observations of TeV blazars are: (i) 
the study of the matter, energy content, and structure of jets from Active 
Galactic Nuclei (AGNs), (ii) the time resolved investigation of particle accel- 
eration processes, (iii) the research on the connection between the accreting 
supermassive black hole system and the jet parameters, and (iv) the measure- 
ment of the total infrared/optical luminosity of the universe. 

In this paper we give an update on the status of the observations (Sect. 2) 
and the four major scientific topics (Sects. 3-6). We conclude with an out- 
look in Sect. 7. We will focus on observations and modeling of TeV blazars 
(blazars detected by TeV instruments). Reviews on observations and models 
of MeV/GeV blazars (detected by EGRET) have been given in (64; 2^; recent 
overviews of TeV astronomy have been given in (jsst iisl ) . 

2 Observations 

TeV 7-ray observations have achieved a state of maturity with good agreement 
in fiuxes and reasonable agreement in energy spectra between different exper- 
iments. The measurements of the steady flux from the Crab Nebula with the 
CAT, HEGRA, and Whipple experiments lie within 18.5% from the common 
mean comparable to the systematic error of about 22% quoted by the groups 
(IstI : 0; Is^ l. The spectral indices agree to within 0.2 from the common mean, 
a difference that is larger than the systematic errors of between 0.04 and 0.06 
estimated by the 3 groups. Systematic error on flux and spectral changes from 
one day to the next are estimated to be about 10% and 0.05, respectively. 

High-quality TeV 7-ray energy spectra have now been measured for Mrk 421, 
Mrk 501, H 1426-1-428, and lES 1959+650. Fig. 1 compares the X-ray and 
TeV 7-ray SEDs of these 4 sources. Compared to the other 3 blazars, the 
low-energy (synchrotron) component of Mrk 421 peaks at lower frequencies. 
Comparison of the high-energy TeV 7-ray energy spectra is not straight for- 
ward as it depends on the uncertain extent of extragalactic absorption (see 
Sect. 6). Spectral variability has been established for the two strongest ones, 
Mrk 421 Hi) and Mrk 501 (011). The significance of measurements of spec- 
tral changes can not be overestimated. Spectral changes (like fiux changes) are 
free of the still substantially uncertain extent of extragalactic 7-ray absorption. 
Modeling of simultaneously taken X-ray and 7-ray data with well measured 
spectral changes in both bands has not yet been performed but has the po- 
tential of breaking basically all model degeneracies. 
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Fig. 1. Simultaneous and 
non-simultaneous X-ray 
and TeV 7-ray energy 
spectra of the 4 TeV 
blazars with measured 
TeV 7-ray energy spec- 
tra. The regions show 
the range of values that 
have been observed with 
BeppoSAX, RXTE and 
Cherenkov Telescopes 
(from 



Several extensive multiwavelength campaigns have been performed. The cor- 
relation of X-ray and TeV 7-ray fluxes had first been seen in Mrk 421 data 
(E3;|63). In the meantime, correlations with high statistical significances have 

I). A more re- 



been observed from Mrk 501 ()45l : l63f l and from Mrk 421 
cent campaign on the TeV blazar lES 1959+650 shows an orphan 7-ray flare 
without an X-ray or optical counterpart (jifih . The orphan flare shows that the 
simplest models (1-zone Synchrotron Self-Compton models with a randomly 
oriented magnetic flelds and isotropic electron distributions) are not able to 
account for the data. 

Simultaneous measurements of blazar energy spectra with imaging Cherenkov 
telescopes and solar array Cherenkov detectors have the potential of produc- 
ing highly interesting results. First lightcurves from STACEE fllil ) and a flrst 
energy spectrum from CELESTE (58) have been reported recently. 

3 Emission Models and Implications for the Jet Structure 
Synchrotron-Compton Models: In Synchrotron Compton (SC) models of 
TeV blazars the radio to X-ray emission is produced as synchrotron radiation 
from a population of non-thermal electrons (and possibly positrons). The same 
electron population emits 7-rays through Inverse Compton processes by elec- 
trons scattering synchrotron photons in so-called Synchrotron Self-Compton 
(SSC) models or with "external" photons that originate outside the jet in so- 
called External Compton (EC) models (see (jc^) and references therein). In the 
case of BL Lac objects, the lack of strong emission lines is commonly taken as 
evidence that ambient photon flelds are not important, and that SSC models 
rather than EC models are more likely to explain the data. 

A method to distinguish between SSC and EC models is to measure the time 
lag Ats_ic between X-ray and 7-ray flux variability (24). In SSC models one 
expects a time-lag of the order of R c^^ with R the radius of the emission 
volume and c the speed of light. The time lag originates from the fact that a 
population of freshly accelerated electrons immediately radiates synchrotron 
emission; in contrast, the Inverse Compton component peaks only after the 
new synchrotron seed photons had time to propagate through the source. No 
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such time lag has yet been observed with high statistical confidence. 

TeV blazar data have been modeled intensively with a variety of codes. The 
more simpler ones assume 1-zone homogeneous emission volumes and use 
phenomenological electron energy spectra to fit the broadband data (e.g. 
fliil : liii fl^) ■ Modehng with time- dependent self-consistently evolved codes 
(e.g. liJE; 44; 3)) assures that the energy spectra are physically realizable 
from initial acceleration energy spectra. Time dependent modeling has to be 
used whenever the fiare duration is comparable or shorter than the longest 
time scale of the microprocesses that modify the electron energy spectra (e.g. 
acceleration energy gains, radiative energy gains and losses). 

Several conclusions can be drawn from the SSC and EC modeling performed 
by various authors: (i) relatively simple 1-component and 2-component SSC 
models are able to explain most multiwavelength data; (ii) even for the nearest 
blazars with z ~ 0.03, extragalactic absorption has to be taken into account 
and substantially modifies the regions in the parameter space that correspond 
to acceptable solutions (e.g. 0- |9|)'): (iii) the data prefer high jet Doppler fac- 
tor of between 50 and 100 (15 It l47r: (iv) in the jet frame the energy density of 
non-thermal particles is much higher than the energy density of the magnetic 
field (the jet plasma is a low cx-plasma) (e.g. (4^ 3)); (v) only a small fraction 
(< 10~^) of the energy transported by the jet is converted into radiation (e.g. 

The last two conclusions constrain the structure of the jet at its base. While 
electromagnetic models of jet formation produce high-cr jets with large electro- 
magnetic field energy to particle energy densities, SSC models indicate low-cr 
jets. The low radiative efficiency indicates that the sub-parsec jet can smoothly 
evolve into the parsec-scale jet, without the need for plasma re-coUimation or 
re-acceleration that would have to follow a more efficient conversion of bulk 
plasma energy into random particle energy and radiation (see also (|2ol)). 
Hadronic Models: In hadronic models the continuum emission is explained 
by hadronic interactions of a highly relativistic baryonic outfiow which sweeps 
up ambient matter fj59l) , by interactions of high-energy protons with gas clouds 
moving across the jet |25|) . or by interactions of Ultra High Energy protons 
with ambient photons (pof). with the jet magnetic field (jlj), or with both (js^)- 
Although attractive since they involve Ultra High Energy Cosmic Ray acceler- 
ation, hadronic models have difficulties to explain the good X-ray/TeV 7-ray 
fiux correlations found for Mrk 421 and Mrk 501. Furthermore, little work has 
yet been done to see if proton models can explain more than single snapshots 
of data. 

4 Particle Acceleration 

Among the many particle acceleration mechanisms that are discussed, Fermi 
acceleration at strong Magnetohydrodynamic (MHD) shocks has received most 
attention. Indeed, the data are largely consistent with an electron acceleration 
energy spectrum dN/dE oc with p between 2 and 2.23, expected for non- 
relativistic and ultra-relativistic shocks, respectively (see ()4lh and references 
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Fig. 2. Correlation of estimated black hole mass and the peak of the Synchrotron 
SED (a) and observed TeV 7-ray flux variability time-scales (b) for the 6 well 
established TeV blazars. Black hole masses from stellar velocity dispersion are shown 
by solid (from (j^) and dashed (from |2^) error bars. Black hole masses from the 
galactic bulge luminosity (from (|23)) are shown by dotted error bars. Vertical error 
bars show the range of observed values. Horizonthal error bars are statistical errors 
in the case of solid and dashed lines, for the dotted estimates no statistical errors 
have ben published and we assumed Alog{M,) = 0.25 (figures from (46)). 

therein). Based on a time-dependent SC code Kirk & Mastichiadis (1999) 
have shown that Fermi energy gains and radiative energy losses can result in 
characteristic hardness-intensity correlations during individual flares (^3)- X- 
ray and recently also TeV 7-ray data do show the predicted signatures: during 
flares the sources seem to go through clockwise and anti-clockwise loops in the 
X-ray or 7-ray hardness-intensity planes (e.g. (0;lai)). However, the same 
source shows a wide range of dffiereut signature7(e.g. B), .ud.cafug that 
the relative length of the characteristic time scales of particle acceleration and 
radiative cooling change from flare to flare. The lack of a prevailing signature 
has cast doubts on whether "cooling loops" and "particle acceleration loops" 
have really been observed. 

Most SC models of TeV blazars invoke a minimum Lorentz factor 7n,in of ac- 
celerated particles on the order of 10^ to account for the data (e.g. (j44f ) and 
references therein). Such high 7inin- values could either result from a yet un- 
known pre-acceleration mechanism, or, from an upstream bulk Lorentz factor 
on the order of 10^. A low-entropy, extremely relativistic upstream plasma 
could be detectable by its Inverse Compton emission if exposed to a suitable 
external radiation fleld. 
5 Black Hole - Jet Connection 

From the theoretical side, the connection between black hole and jet proper- 
ties is far from being understood. Purely electromagnetic |l^ l49l l and GR- 
electromagnetic models (0) of jet formation face the well known a-problem, 
namely that they predict Poynting flux dominated energy transport while 
SC models indicate strongly particle energy dominated iets. Hydrodynamic, 
MHD, and GR-MHD scenarios of jet formation (e.g. (0 EJ) face the diffi- 
culty that they predict jets with bulk Lorentz factors F < 10 while SC models 
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favor r > 10. 

From the observational side, first studies of the black hole-jet connection are 
just becoming possible based on high-quality black hole mass estimates from 
galactic stellar velocity dispersion measurements ([iii H). The studies per- 
formed so far do not indicate any correlation of the black hole mass and the 
peak location or peak luminosity of the low-energy or high-energy SEDs, the 
time scale of flux variability (see Fig. 2), or the flare duty cycle. The cor- 
relations may be masked by variation of parameters like jet viewing angle, 
jet mag netic field, or intensity of the amHent photon field i M). Alterna- 
tively, the black hole mass estimators may not work in the case of blazars, as 
discussed in (jsl). 

6 Measurement of the Cosmic Infrared Background 

Extragalactic absorption of TeV 7-rays by the Cosmic Infrared Background 
(CIB) and the Cosmic Optical Background (COB) in 7Tev + 7cib/cob 
e"*" e~ pairproduction processes ((5^ ,31,; IgA lioh substantially complicates 
the study of the astrophysics of blazar jets. However, the absorption effect 
allows us to use TeV 7-ray observations to measure the energy spectrum of 
the CIB and the COB. These backgrounds are extremely interesting owing 
to their close relation to the total electromagnetic luminosity of the universe 
since decoupling of matter and radiation ~300,000 years after the big bang. 

The TeV energy spectra of extragalactic sources carry the imprint of extra- 
galactic absorption in the form of "characteristic high-energy cutoffs" in the 
30 GeV to 10 TeV energy range. The measurements of these cutoffs makes 
it possible to constrain the energy spectrum of the infrared to optical back- 
ground radiation. These constraints are very valuable, as direct measurements 
are plagued by large systematic errors, owing to strong foreground emission 
from our solar system and the Milky Way (e.g. (H^). While the spectral mea- 
surement of the nearby blazars Mrk 421 and Mrk 501 {z ^ 0.03) have been 
used to set upper limits on the CIB intensity in the 1-20 /im wavelength range 
( 111), the recent measurement of the energy spectrum from the more distant 
blazar H 1426+428 {z = 0.129) has given the first demonstration of the full 
potential of TeV 7-ray observations to decide between different background 
models 

Further progress may best be accomplished with a two-pronged approach: (i) 
Detailed source modeling makes it possible to disentangle, at least partially, 
source inherent high-energy cutoffs from those caused by extragalactic extinc- 
tion and thus to improve on the constraints from each single source (23; .44) : 
(ii) since the spectral cutoffs from extragalactic absorption are expected to 
occur at energies which depend only on source redshift, a statistical analysis 
of the constraints from many sources can be used to test and strengthen the 
constraints from individual sources. 

7 Outlook 

Over the last several years, the field of TeV blazar studies has benefited greatly 
from detections of new sources and from very intensive multiwavelength cam- 
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paigns. Further progress may derive from extending the time scales of the 
broadband observations. High-accuracy measurement of Gamma-Ray Burst 
(GRB) afterglows is boosting the theoretical progress in the case of GRBs. 
In the same way, detection of increased flux levels in the radio to UV bands 
following major blazar flaring phases may prove crucial for pinning down the 
jet properties. For this purpose regular broadband monitoring is required over 
time scales longer than the typical 10 days of current multiwavelength cam- 
paigns. As the quality of the data improves, it now becomes apparent that 
simple one- zone SC model fit snapshots of the data, but are not able to give us 
deep insights into the jet structure and the nature of the particle acceleration 
processes. Most present theoretical work focuses on the emission mechanism 
and largely neglects the flare origin and the development of the plasma that 
supports the non-thermal particles during the flare. A notable exception is 
the recent work concerning the internal shock model (0;0t). In future work, 
combining 3-dimensional MHD simulations of flares with kinetic simulations 
of the non-thermal particle population would make it possible to explore more 
complex observational signatures than possible with present codes. 

In view of the next-generation Cherenkov telescopes CANGAROO III, H.E.S.S., 
MAGIC, and VERITAS presently under construction, the upcoming launches 
of the X-ray observatories SWIFT and ASTR0-E2, and the 7-ray observato- 
ries AGILE and GLAST to be launched in 2005 and 2006, respectively, we 
can expect a wealth of exciting blazar results in the next 5 years. 
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